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We perform a comparative ab-initio study of Mn-doped Germanium semiconductor using 
the Perdew-Burke-Ernzerliof (PBE) exchange-correlation functional, DFT-|-i7 and Heyd-Scuseria- 
Ernzerhof hybrid functional (HSE). We show that the HSE functional is able to correctly account 
for the relevant ground state properties of the host matrix as well as of Mn-doped semiconductor. 
Although the DFT+L'' and the HSE description are very similar, some differences still remain. In 
particular, the half- met allicity is lost using DFT+U when a suitable U value, tuned to recover the 
photoemission spectra, is employed. For comparison, we also discuss the case of Mn in Silicon. 
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INTRODUCTION 



Dilute magnetic semiconductors (DMSs) are still a 
topic of great current interest. The theoretical de- 
scription of the interaction of transition metal doped 
semiconductors is challenging since localized states in- 
teract significantly with delocalized states. Density func- 
tional theory (DFT) in the local density or generalized 
gradient approximation (LDA or GGA) for the exchange- 
correlation energy is not able to properly describe the 
non-locality of the screened exchange interaction and, 
furthermore, possesses a sizeable self- interaction error. [5j 
These limitations are particularly severe in the case of lo- 
calized orbitals, e.g. Mn-Sd states, which are described as 
too shallow in energy resulting in a large hybridization 
with anion p-states. As a result, the Mn-Sd states are 
over-delocalized. The situation is particularly serious in 
the case of small band-gap semiconductors (such as Ge) 
which are described as metals in LDA/GGA, thus pro- 
ducing an overestimated hybridization among the valence 
and conduction Ge states and Mn-d states. The physics 
of localized d states can be partially described using a 
'DFT+U formalism, which introduces a local correction 
U to recover the proper position of the Mn-d states Q. 
However, in Ge, the accurate electronic properties are 
not completely recovered: e.g. the half-metallicity of the 
compound is lost within DFT+U scheme. 

Very recently, hybrid Hartree-Fock density function- 
al, which mix a fraction of the exact Fock exchange 
with the DFT exchange, have been widely applied to 
extended solid state systems. 0, 0-23 1 In this paper, we 
mainly focus on Mn-doping in bulk Ge by performing 
hybrid-density functional theory calculations. We will 
show that the HSE functional gives a satisfactory de- 
scription of the structural, electronic and magnetic prop- 
erties of Ge-based DMS, consistent with experimental 
data. Furthermore, for few selected properties, and, for 
sake of comparison, we will also include some results of 



Mn-doped Silicon. 



COMPUTATIONAL DETAILS 

The calculations were performed within the projector 
augmented-wave (PAW) methodj24i] using the Perdew- 
Burke-Ernzerhof (PBE) generalized gradient approxi- 
mation (GGA)'25'| and Heyd-Scuseria-Ernzerhof (HSE) 
hybrid functional [2^, \2l\ . recently implemented in the 
VASP code.jii, also used the DFT+U method 

for Mn within Dudarev's approach Is^l fixing J7 to 6 eV 
and J to 1 eV (PBE was used for the DFT part). The 
kinetic energy cutoff used for the orbitals was set to 300 
eV. Monkhorst-Pack fc-point grids of 10 x 10 x 10 and 
6x6x6 were used to sample the Brillouin zone of the 
Ge-bulk and of the 64-atom unit cell, respectively. All the 
atomic internal positions were relaxed. In the following, 
we will focus on the Ge bulk system, single (substitu- 
tional and interstitial) and double substitutional (dimer) 
Mn impurities in a 64-atom Germanium cell. For Silicon, 
we will consider the bulk case and the Mn substitutional 
impurity in a 64-atom unit cell. 



BULK GE 

For bulk Ge, the calculated equilibrium properties are 
in good agreement with experiments: the HSE lattice 
constant is 5.703 (5.792) A using HSE (PBEVwithin 0.7 
(2.3) % of the experimental value of 5.660 A;[3l| the HSE 
bulk modulus (731 kbar) improves over the PBE value 
(571 kbar) when compared to experiment (768 kbarfsJ]). 
Furthermore, we remark that within HSE the energy 
gap is properly described to be indirect (0.63 eV includ- 
ing SOC[33^ compared to the experimental value of 0.74 
eV|34]). A similar result was obtained within HSE for 
Sil33|: the HSE lattice constant is 5.444 A using HSE[3l| 
compared to the experimental one of 5.430 A. [^] Thecal- 



is 1.12 eV[33[ while the exper- 



culated indirect energy ga 
imental one is 1.17 eV.IS^ 



It is interesting to note that, 
for Silicon, self-interaction schemes which are often used 
to improve the electronic structure description,!^ [s^ do 
not open the gap.fsl] This is important for our present 
study, since a faithful description of the equilibrium prop- 
erties of the host semiconductor is at the basis of an ap- 
propriate description of the doped system. 



MN IMPURITIES IN GE 

In TabUl we summarize our main results, i.e. forma- 
tion energyls^ (AH), Mn-Ge bond length (divin-Go) and 
the Mn magnetic moment (//) at their respective theoret- 
ical lattice constants for the considered Mn-doping cases. 
The formation energies are evaluated with respect to the 
calculated Ge and Mn equilibrium bulk phases (diamond 
Ge and AFM-fcc Mn). Of course, while Ge-rich growth 
conditions can be safely assumed to fix the Ge chemical 
potential to its bulk value, the same is not true for Mn so 
that the formation energy is a function of the Mn chem- 
ical potential ^lun- In Table U we report the value for 
Mn-rich conditions fixing to the corresponding bulk 
value (a-Mn). 

We note that the experimental evidence of local Ge- 
lattice dilation upon Mn-doping is correctly described 
using HSE yielding a Mn-Ge distance 2 % larger than 
the ideal Ge-Ge bond-length, while PBE gives a local 
contraction of —2% and DFT-I-J7 finds a smaller local di- 
lation of 0.4 (0.8)% when using the theoretical PBE (HSE 
or experimental) lattice constant. The most recent 

extended x-ray absorption fine structure results [ill - fisj 
yield a Mn-Ge coordination distance of 2.50-2.51± 0.03 
A for the samples obtained at low temperature and which 
are thought to be best candidates for Mn occupation on 
substitutional sites [i^. Clearly, these results match the 
HSE result and the DFT+U as well. Similar results for 
bond-lengths contraction/dilation within different DFT 
schemes were reported also for III-V based DMS[44]. 

A simple molecular orbital description, as sketched in 
Fig. [U can be useful in order to describe the interaction 
of Mn in the tetrahedral Ge ligfmd field, as also done pre- 
viously for similar compounds [45j. We recall that in dia- 
mond like semiconductors, the sp valence states arrange 
to form sp^ hybrid orbitals, each of them filled with a 
bonding electron pair. If one Ge atom is removed creat- 
ing a Ge vacancy, 4 sp^ hybrids point towards the vacant 
Ge atom, each filled with one electron (dangling bonds). 
The Ge vacancy is now replaced by a Mn atom. Due to 
the local tetrahedral symmetry, the Mn d-states are split 
into 3-fold degenerate t2g and 2-fold degenerate e^-like 
states, further splitted by the local exchange field (see 
Fig. [l] right part). From linear combinations of the four 
sp'^ Ge dangling bonds pointing towards the transition- 
metal impurity, an s-like ai-f.i orbital and three p-like 
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FIG. 1: (Color on line) Molecular energy diagram of the Mn- 
d states (right) interacting with the Ge-sp^ hybrid orbitals 
(left). The Mn-s states are not shown for clarity (see text), b, 
ab and nb subscripts label bonding, anti-bonding and non- 
bonding orbitals. Arrows denote up/down electrons while 
circles indicate holes. The central panel (Mugg) shows p — d 
hybridization for substitutional Mn at Ge site. 



t2gt,i orbitals are formed. The host ai-f,^ orbital and 
the transition-metal 4s-f-4 states form a doubly occupied 
bonding state deep in the semiconductor valence band, 
and an empty antibonding state high in the conduction 
band (not included in Fig. [!}. For the majority compo- 
nent, the yhi-t^g orbitals are lower in energy than the 
Ge-t2g sp^ hybrid states. They interact giving rise to 
3 bonding states (3xMn-i2g)6 (see Fig. [T]) and 3 anti- 
bonding states (3xGe-sp'^)ab. The Mn-Cg states do not 
hybridize because they are non-bonding in a tetrahedral 
ligand field. For the minority component, the Ge-t2g 
sp^ orbitals are lower in energy than the Mn-t2g states. 
Upon interaction, they give rise to 3 bonding orbitals 
(3xGe-sp'^)b and 3 anti-bonding orbitals (3xMn-i2g)afc- 
The complex is characterized by a total of 11 electrons 
(4 from the nearest Ge atoms and 7 from the Mn im- 
purity atom). Disregarding the two electrons occupying 
the lowest oi symmetry-like state, one needs to fill the 
orbitals with nine electrons as shown in Fig. [T] Clearly 
the Mn impurity is in a high spin state with 5 d-electrons 
in the majority channel (2 electrons in the Cg and 3 in the 
t2g-like states) and zero d-electrons in the minority chan- 
nels. However, while the minority valence Ge-sp^ states 
are fully occupied, the majority states accommodate two 
holes. This simple molecular picture suggests that: i) the 
compound is half- metallic, with the Fermi level falling 
within the Ge majority valence band; ii) the total spin 
moment of the complex, i.e n-]- — n^, is 3 /is; iii) the local 
Mn d spin moment is 5 /is partially compensated by the 
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holes in the sp^ states; iv) the induced spin moment on 



the 4 nearest Ge atoms, i.e. n^°~*^ — n^'^~''^ should 
be sizeable and opposite to the s pin on the Mn atom. 
In line with previous calculations 46 - 4^ , the calculated 
results (see Tab. HI confirm this picture finding a total 
spin moment of exactly 3 /is in the unit cell, 4.1 fj,B at 
the Mn atom, and —0.11 at the nearest Ge-atoms. 
Furthermore the sizeable induced moments on Ge atoms 
suggest that the holes are rather delocalized. 

The local angular momentum decomposed density of 
states (DOS) shown in Fig. [5] is consistent with this or- 
bital interaction diagram. The top panel shows the HSE 
results, whereas the bottom panel reports DFT+C/ re- 
sults for U=6 eV. In the inset, we show the relation be- 
tween the center of mass of the Mn-d majority states, 
(ed), and the U value. The horizontal line indicates the 
(cd) value, which matches the HSE result {U=6 eV). 
Fig.[2]clearly confirms the interpretation discussed above. 
In particular, integration of the majority total density 
of states from the Fermi level up to the end of the Ge- 
valence band exactly sums up to 2 electrons: these are 
the two holes required to fill the Ge valence band. Mn- 
substitution into the Ge host matrix does not produce a 
Jahn- Teller ion, but rather a Mn^+ ionic state with a 
configuration and two spin-polarized holes. Let us now 
compare the HSE and BFT+U DOS. Due to the choice 
of the U value, the Mn-d states have the same energy 
within HSE and DFT+U but the hybridization between 
Ge-sp and Mn t2g and Cg states (this latter symmetry- 
allowed away from F) is underestimated within DFT+U 
compared to HSE. While the on-site U mainly localizes 
the Mn d states, HSE also acts through the screened 
exchange on Ge-p states, lowering their energy position 
and leading to a larger Mn-d Ge-p hybridization. As a 
matter of fact, the larger hybridization in HSE compared 
to DFT+U can be recognized just below —4 eV, where 
a peak in the eg character (shaded region in Fig. [2|) is 
completely absent in the present DFT+U description. 
We note that previous LDA+U calculations^ with U =4 
eV, although reproducing the peak at —4 eV character- 
istic of the Mn-Ge bond, Q gave a quite different density 
of states for both the t2g and eg states, which is due to 
the strong dependence of the localized d-states descrip- 
tion on the U parameter. Finally, as found in Ref. |a, 
the half-metallic character of the compound is destroyed 
within DFT+U: the energy position of the Mn-t2g Ge-sp 
bonding minority states, whose energy position is mainly 
determined by the atomic Ge-sp levels, is raised towards 
higher energies causing an incomplete filling of the mi- 
nority valence band. 



In a previous study, [47[ it was shown that the half- 
metallicity is favored in Mn doped Ge while in Silicon 
matrix it is lost. In Fig. [3l we show the HSE DOS for 
substitutional Mnsi (top panel) and DFT+U, with U =6 
eV, the same used for Muqc (bottom panel). The cor- 



TABLE I: Formation energy AH, Mn-X distances dMn-x and 
magnetic moments ^ for Mn-doped Ge for various structures. 
Distances in parentheses specify the ideal Ge-Ge bond-length 
of the host (lines dMn-x)- Local Mn magnetic moment as well 
as total magnetic moment (in parentheses) are specified. 





PBE 


DFT+U 


HSE 


Mn substitutional site 








AH (eV/Mn) 


1.5 




0.9 


dMn-Ge (A) 


2.46 (2.51) 


2.53 (2.51) 


2.52 (2.47) 




3.3 (3.1) 


4.1 (3.4) 


4.1 (3.0) 


Mn interstitial site 








AH (eV/Mn) 


2.1 




1.8 


dMn-Gc (A) 


2.57 (2.51) 


2.63 (2.51) 


2.58 (2.47) 


M (ms) 


3.4 (4.0) 


4.2 (4.8) 


3.8 (4.1) 


Mn-Mn dimer 








AHafm (eV/Mn-pair) 


2.9 




1.5 


AHfm — AHafm 


0.77 


0.22 


0.15 



d, 



(A) 

"Mn-Mn(AFM) (A) 



Mn-Mn (FM) 



2.55 (2.51) 2.97 (2.47) 2.99 (2.47) 
1.95 (2.51) 2.71 (2.47) 2.83 (2.47) 



rection for the self-interaction error has a larger effect on 
Si-sp'^ states, since they are quite localized. Therefore, 
they are pushed down in energy. According to the or- 
bital energy diagram shown in Fig. [1] also the minority 
bonding (Si-sp'^)h are shifted down in energy, favoring the 
half-metallicity. Obviously, the Mn-d states are also cor- 
rected for the self-interaction error, and they are pushed 
down in energy as well On the other hand, DFT+U cor- 
rects only the Mn-d states, but not the Si-sp"^ states. This 
gives a near-half metallic structure and an underestima- 
tion of the hybridization of Mn-Si states compared to the 
HSE description. 



The single interstitial impurity 

This defect in a Germanium matrix possesses twice the 
formation energy as the substitutional impurity (HSE), 
hence it is unlikely to form. Here we only note that a 
tendency to a local expansion around Mn is found for 
all functionals. Interestingly, the local magnetic moment 
is larger for DFT-I-J7 than for HSE, suggesting sizeable 
differences in the interaction of the impurity with the 
local environment. 



Mn-dimers 

Double Mn substitutions on two nearest-neighbouring 
Ge sites are of particular interest, since they are in- 
ferred to occur at experimental growth conditions 41|.|48| 
leading to nucleation of Mn-precipitates. In addition. 



they might be also detrimental for the magnetic order- 
ing, since dimers show antiferromagnetic (AFM) coupling 
with no net spin moment. Unfortunately, GGA-based 
calculations are not entirely conclusive, since the dimer 
configuration becomes stable at a too small bond-length 
(1.95 A) not compatible with the Mn ionic radius (~ 1.1- 
1.3 A) or bond distances in Mn-Ge compounds. Thus, 
published results [H, |4§| often refer to the ideal unre- 
laxed structure which, of course, strongly overestimates 
the heat of formation of the dimer. From the results re- 
ported in Tab. U assuming thermodynamic equilibrium, 
we can comment on the relative concentration of single 
substitutional sites and dimers. At thermodynamic equi- 
librium, the concentration c of a defect with riMn Mn 
atoms is roughly proportional to e^^'^^^"^!"^'^"")/'^^-'" 
where AH is the formation energy, the Boltzman con- 
stant, and T the temperature. Supposing A/j,Mn ~ eV 
(thermodynamic equilibrium with a— Mn), the probabil- 
ity of finding a substitutional Mn or dimer is q-O-^/^^bT 
and ^ respectively, i.e. monomers are more 

likely to form than dimers. More generally, for a spe- 
cific A/iMn the probabilities for single substitutions and 
dimers are, 



02 

o 

G 




,(--0.9+A/iM„)/feBT 



and 



,(-1.5+2MAM„)/fcBT 



Therefore the dimer concentration is larger than the 
monomer concentration only for A/^Mn > 0.6 eV, i.e. 
at extremely Mn-rich conditions, where a— Mn precipi- 
tates are anyway already preferred over the formation of 
monomers (or dimers). Although kinetic effects might 
well hinder the nucleation of larger precipitates, our cal- 
culated thermodynamics suggests a rather low dimer con- 
centration. It is important to note, that the thermody- 
namic arguments alone, presented here, are not enough 
to fully discuss the relative probability of occurrence of 
the monomers with respects to dimers, as these systems 
are usually grown out of the thermodynamic equilibrium 
and kinetic effects may have an important role. 

Finally, the stabilization energy of AFM over FM cou- 
pling of nearby impurities is much lower within HSE than 
GGA. Furthermore, for HSE (and BFT+U), the calcu- 
lated Mn-Mn distance for both FM and AFM magnetic 
alignment, are in line with experiments — always report- 
ing local lattice dilation — as well as in line with the 
Mn-Mn distances in the FM Mn5Ge3 compound (vary- 
ing between 2.52 and 3.06 A) [SQtS^]. 



CONCLUSIONS 

In summary, we have performed a comparative study 
of substitutional Mnce by using FEE, PBE-f [/ and HSE 
functional. The main focus is on the differences aris- 
ing from three different treatments of the exchange- 
correlation term, namely the FEE, DFT+U, and HSE. 



FIG. 2: Density of states projected on the Mn impurity site 
(top panel) in symmetry resolved angular momentum com- 
ponents: t2g (dashed line) and eg (shadow) states for Mn 
substitutional impurity. The density of states projected on 
the I — 1 component of the 4 Ge (bottom panel) coordinated 
with the Mn impurity is also shown (solid line). The inset 
shows the Mn-d center of mass as a function of the U value. 
The horizontal line indicates the value found within HSE. 



As well known, the PEE treatment can not describe sat- 
isfactorily the ground state properties of the Mn in the 
host semiconductor matrix. Including the U correction 
at DFT level improves the description. However, some 
differences still remain when compared to HSE. For ex- 
ample, the HSE Mn-d peak position is found at ^ —5 eV 
with respect to the Fermi energy, i.e. in the same en- 
ergy region as observed in photoemission experiments. 
When using the DFT+U method and fixing the U pa- 
rameter in order to recover the experimental d-peak posi- 
tion, the hybridization with Ge-p states around —4 eV is 
underestimated compared to HSE and the experimental 
photo-emission peak.Q Furthermore for the same U, the 
half-metallic character is not predicted by DFT-I-J7. The 
fact that HSE accurately describes the host semiconduc- 
tor and, at the same time, the interaction of localized 
Mn states with the host valence states makes this func- 
tional a valuable approach for studying transition metal 
defects in semiconductors. It is also true that the HSE 
calculations are usually quite more computationally de- 
manding with respect to F)FT+U. Therefore, whenever 
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E-Ej(eV) 



FIG. 3: (Color on line) HSE Density of states projected on 
the Mn impurity site (top panel) for Mn in Si. Labels as in 
Fig.[l DFT+f/, (7=6 eV (bottom panel). 



a compromise between accuracy and computation effort 
is required in the calculations, a preliminary HSE study 
may be useful for choosing an appropriate U value, which 
is often not accessible from experiments. 

This work was supported by the Austrian Fonds zur 
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computing grant at CINECA-HPC center. 
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